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By FRANK LOXLEY GKIFFIN, Williams College. 



§ 1. Circles as Central Orbits. 

1. The Problem. Given any central force which is everywhere finite 
and real, and which is a function of the distance alone, f=f(r), it is known 
that all circles with centers at the center of force are possible orbits. For 
if a particle were to move in such a circle, of radius a, the centrifugal force 
v 2 /a would have to balance exactly the attraction f(a) ; and, conversely, if 
the particle were started at a distance a from the center of force, in a direc- 
tion perpendicular to the radius vector, and with the velocity defined by 
fl 8 — af(a)> the circular orbit would result. 

But in order that a circle whose center is elsewhere than at the cen- 
ter of force, be described as a central orbit, the force must vary according 
to a special law, which is obtainable when the polar equation of the given 
orbit has been written. The object of this paper is (1) to derive such a law 
of force, (2) to ascertain what other circular orbits are admitted by it and 
what properties have they in common, (3) to learn what other laws admit 
families of orbits having some of those properties, and (4) to point out the 
families in question in the case of two simple laws. 

2. Laws Admitting Circular Orbits. Let the center of force be 
selected as the pole, and let the polar angle o be measured from the 
line through the pole and the center of the circle; then the equation of the 
orbit is 

(1) 2dcos 0=l/i*- (a 2 -d*)u, 



where u is the reciprocal of the radius vector, a is the radius of the circle, 
and d is the distance between the pole and center. [It is necessary to take 



•Read before the American Mathematical Society, September 10, 1908. 
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d<a, since a closed central orbit not containing the center of force is 
impossible. *] 

The required central force must vary along the orbit according to the 
law: f=h s u i (u+d ! u/dO i ), where his the constant of areas. The differ- 
entiation of (1) with respect to o t followed by the elimination of o from the 
two equations gives 

(2) W = \l/u i -(.a*-d i )Yu* + \l/u i + {a i -d*)Y i {du/do)\ 
whence is obtained by adding 4(a s — d 2 ) to both members: 

(3) u s + (du/do)^=Aa 2 u i [(a 2 -d i )u'+l']- 2 . 

By differentiating again, or by noting that the curvature, given by 
u 3 (u+d'u/d o s ) [u* + (du/d o) *]-', equals 1/a at all points of the circle, the 
necessary force is easily obtained as a function of u alone: 



(4) f= 



(w 2 w 2 +l) s (w 2 + r 2 ) 3 ' 

where m 2 =8a 8 /?. 2 and n 2 =a 2 — d % . Some obvious properties of this force 
are: (a) it vanishes both "at infinity" and at the center of force; (6) it is 
everywhere real; and (c) it is everywhere finite and continuous. If 
the values of a, d, and h had been different in the given orbit, the values of 
m and n would usually be different; but the properties mentioned would 
persist. The type of law where the force is given by (4) for particular 
values of m and n will be called Type I hereafter. 

§ 2. Other Circular Orbits for Laws op Type I. 

3. A Family of Circles. Consider any particular law of Type I. 
Because of its properties, just enumerated, this law admits the infinitude of 
circular orbits whose centers are at the pole. But, leaving aside this some- 
what trivial family, the question arises as to the possibility of further circu- 
lar orbits besides the original one. Evidently in all such orbits which may 
be admitted by the given law, a*h 2 and a 2 — d* must be invariant, for 
if either m or n changes value, the law changes. 

The constancy of the second expression restricts the set of circles ad- 
missible as orbits; that of the first imposes a restriction upon the initial con- 
ditions. The totality of possible circles is represented by the equation 

(5) 2Dcos(o-K)=l/u-n { u, 

*Cf . E. J. Routh, Dynamics of a Particle, pp. 292-3. 
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obtained from (2) by replacing the original d and a by D and A the corres- 
ponding constants for any other circle, imposing the restriction A 2 — D 2 = n 2 , 
and introducing K as the longitude of the center of the circle. There 
is a double infinitude of these circles, obtained by varying D and K. 

Moreover, as will now be shown, there exist real initial conditions for 
which any one of these circles is described as a central orbit. * 

4. Circular Orbits through a Point. Theorem I. Through any 

point (U[, o^ there passes in each direction one and only one circle of the 

family (5). For, if the given direction makes with the radius vector to that 

point an angle <f> lt then, since tan <P=r do/dr=—u do/du, equation 

(3) gives for any circle of (5) on replacing a by A and a 2 — d 2 by n 2 : 

(6) u?<XBC % 4^=lA % u?[n*u?+\]-*. 

This equation determines a single real positive value of A, which is never 
less than n, since cscVi = 1 and n 2 u 2 + 1 = 2nu\ : 

(7) A 5 =--cscV,.(% 2 % 1 2 +l) 8 /4w 1 8 ±n\ 

Hence D 2 =A* — n 2 gives a real value for D. The substitution of the latter 
value, with the coordinates of (u u °i) in (5) determines if uniquely (except 
for additive multiples of 2 *■) and real. 

Theorem II. There is a unique velocity with which a particle project- 
ed at any point in any direction will describe a circle. For, by the preceding 
theorem, A is determined uniquely; and hence HA—m determines H 
uniquely. And since H=v l r^sm'I'i in any central orbit, there is but one 
possible value for v t . Now to describe the required circle with the 
determined value of H, the particle would have to be attracted by a force 
whose law is precisely (4) ; hence, conversely, projection in the given direc- 
tion with the velocity V\ will result, under the law (4), in a circular orbit. 

5. Properties of the Family of Orbits. From the constancy of 
AH it follows that the constant of areas required in any circular orbit de- 
creases as the radius of the circles increases. From the constancy of 
A 2 -D 2 it follows that the radius of the circle must always increase with 
the distance of its center from the pole. There is no upper limit for A or 
D, but there is a lower limit (not zero) for A: viz., the smallest possible 
circular orbit is that one whose center lies at the center of force. But the 
two invariants have more exact interpretations, which are both simple and 
interesting. 

"In a paper "On the Law of Gravitation in the Binary Systems," [American Journal of Mathematics, vol. 31, 
pp. 62-85], the writer discussed the family of conies by a law of which Type I is a limiting case. The present fam- 
ily of circles, however, is not included among those conies; for each law there considered admitted only one circular 
orbit, and the center was at the pole. Certain questions raised and answered below would not have arisen in con- 
nection with the family of conies. 
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Property 1. The curvature at either apse is proportional to the con- 
stant of areas required for the orbit For the curvature is 1/A which equals 
H/m. 

Property II. The product of the two apsidal distances is constant. For 
these distances are A — D and A+D whose product equals n 2 . 

Property III. All circular orbits whose apsidal radii vectores are col- 
linear pass through two fixed points on the perpendicular drawn to the line 
of centers at the center of force. For the half -chord intercepted on this per- 
pendicular, being a mean proportional between A-D and A+D, is constant 
and equal to n 

Another (and more useful) statement of this property is that the vec- 
torial angle from either apse to the radius vector of length n is constant and 
equal to */2. 

§ 3. Other Laws Admitting Properties I— III. 

6. Admitting Property I. Two of these three properties were 
found also in the family of conies mentioned above; and it is natural to in- 
quire whether there are still other laws of force whose trajectories will pos- 
sess such properties. Attention will be confined to those cases where the 
force is a continuous function of the distance. 

Let the force be denoted by u 2 P(u), and let F(u)=2fP(u)du. 
[Since the continuity of P(u) follows from the hypothesis that /is a contin- 
uous function of r, P(u) is integrable; thus F{u) exists and is continuous.] 

A first observation is that any such force, whether attractive 
or repulsive, admits families of orbits having pericenters at arbitrary 
distances; and, if everywhere attractive, admits also families having apo- 
centers at arbitrary distances. For a particle may be so projected at any 
distance 1/P that the initial direction is perpendicular to the initial radius 
vector, so that du/d 0=0 for u—P\ and the initial velocity may be such as to 
give h any desired value. Moreover the standard differential equation of 
central orbits, 

(8) f=h 2 u*(u+d 2 u/do'>)=u*P(u), 

shows that d*u/d o* > at u— P according as h 2 % P(P)/P. Hence whatever 
the sign of P(/5), h may be chosen so large that d i u/do t <Q, making u a 
maximum at the initial point. And, if P(P) >0 [force attractive], it is also 
possible to choose h so small that u is a minimum at u=P. 

Secondly, since the expression for the curvature [see No. 2] becomes 
at an apse u+d^u/d o 2 , or simply P(u)/h 2 , Property I would require: 
h=iP{P)/h*, or 

(9) h*=KP(P), 
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where * is a constant factor of proportionality. Thus for any value of *, h 
is defined as a function of P so as to insure Property I. It remains to be 
shown simply that -* may have such values that h, defined by (9), will meet 
the condition above for the existence of a pericenter or apocenter at u=P. 
Now to have h*>P(P)/P is to have [*P(P)] 2 > [P (P)/P] \ or **>P(P)/P 3 . 
But for any interval of P, 0<A<P = B, P(P) being continuous has an upper 
bound M. Hence, if * 2 be chosen greater than M/A a , the condition is satis- 
fied. And if P(P) >0 always, the reverse inequality may be treated in like 
manner. 

Therefore, every central force which is a continuous function of the dis- 
tance admits families of orbits having pericenter s and possessing Property I; 
and, if everywhere attractive, families of orbits with apocenters, having that 
property. 

Forces admitting orbits with both a pericenter and an apocenter are 
somewhat special as will appear in the next paragraph. For the orbits to 
have Property I at their pericenters only, any continuous force whose tra- 
jectories have two apses is admissible, the corresponding apocenters being 
determined as follows: One integration of (8) gives 

(10) h i (u i + (du/dO)*)=c+F(u), 

where c is an arbitrary constant. If P and « denote, respectively, the peri- 
central and apocentral values of u, then 

(11) h 2 {P 2 -« 2 )=F{P)-F{*). 

Now if P be allowed to vary, and h be determined by (9), then (11) will de- 
termine « in those cases where the orbit has two apses. If, however, Prop- 
erty I is to hold for both apses, there is a further restriction upon the force, 
given by ■* P(P) =/* P(«), a* being the factor of proportionality for apocenters. 
A case of this sort is treated later [in No. 9]. 

7. Admitting Property II. If an orbit has two apsidal distances, 
l/« and VP, Property II requires simply that ap—k*, a constant As 
remarked, however, there is a second category of forces none of whose tra- 
jectories have two apsidal distances; an example is the case where the force 
varies inversely as the cube of the distance. Still other laws admit trajec- 
tories consisting of two branches separated by a region of imaginary motion, 
—one branch having only a pericenter and extending to infinity, the other 
having only an apocenter and reaching to the center of force. In such a 
case which branch is actually described in the motion of the particle 
depends upon the initial distance of the particle; and the path really has but 
one apse. It is easy to find conditions upon the force, which are necessary 
and sufficient to place it in any of the three categories. 
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Theorem I. In order that a force, f, shall admit orfyits extending from 
one apse to another, it is necessary and sufficient that f/u' 6 be a decreasing 
function of u throughout some interval of u, say [a, 6] . For if there is an 
apocenter at %=« and a pericenter at u=? [£><*], the necessary constant of 
areas is given by (11); and moreover it is necessary that h s d' 2 u/do 2 , 
or — h 2 u+P(u), be positive for w=« and negative for u—P; thus 

(12) P(a)/«>V>P(/»)//S. 

If P(u)/u is nowhere a decreasing function of u, it is clear that the inequal- 
ity (12) is impossible; thus the hypothesis is necessary. 

It is also sufficient; for let a =a and /?=&. Then the use of the gener- 
alized theorem of mean value in (11) gives the constant of areas for which 
du/d o vanishes at w=« and at u=P: 

as) h .jm=m = p(2 ll 

where *<Ui<P, the theorem being valid here, since F(u) and u 2 together 
with their derivatives are continuous functions in [«, /?], and the latter de- 
dominator does not vanish in the interval. Now by hypothesis, 

(14) m_ > PM > m_. 

v ' a U, P 

so that (12) is satisfied, and there is an apocenter where u—a and a peri- 
center where u=P. Whether the orbit extends from u =« to u—P or wheth- 
er there are still other apses between these values of u (for the same values 
of h) does not affect the theorem. For consider the next apsidal value above 
u= a , say u=". Then since h 2 =P(u 2 )/u s where <*<«*<«'; and since this 
shows that h 2 >P(<*')/ a ', the apse at «=«' is a pericenter. A similar argu- 
ment shows that the apse nearest to u=P, say u=P<P, would have to be 
an apocenter. In any event the same branch has both a pericenter and an 
apocenter. 

Corollary 1. If f/u 9 decreases everywhere, then any two distances 
can be the apsidal distances in an orbit. For the constant of areas required 
for the vanishing of du/d o would be such as to give a pericenter at the les- 
ser distance and an apocenter at the greater. 

Theorem II. In order that a force, f, admit trajectories consisting of 
two branches, one having an apocenter only, and one a pericenter only, it is 
sufficient that f/u 3 be everywhere an increasing function ofu, and necessary 
thatf/u % bean increasing function in some interval [a, &]. For the exis- 
tence of a pericenter at u— P and an apocenter at u=" [ a >P] requires the 
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inequality (12), which can be satisfied only if flu 3 is an increasing function 
in some interval between P and «. Further, since (14) is satisfied iif/u 3 is 
everywhere an increasing function, the apses at u=* and u=P must be, re- 
spectively, an apocenter and a pericenter; and there are no further apses by 
Theorem I, the necessary condition not being fulfilled. 

Corollary 2. If flu 3 is constant [the case of the inverse cube of the 
distance] , only one apse is possible in a trajectory. For the conditions nec- 
essary for two apses either in the same branch or in different branches are 
not satisfied. 

Corollary 3. There is no apsidal distance between l/« and IIP in the 
orbits of Theorem I. For between the apocentral value u—P' and the peri- 
central value w=«', flu 3 would have to be somewhere an increasing function. 

Theorem HI. Any law admitting trajectories with two apsidal dis- 
tances admits families having Property II. 

Case I. When there is an interval [a, 6] throughout which f/u 3 de- 
creases, any two apsidal values « and P [a^«</? < b] determine a value of 
h for which the orbit lies in the region u=a to u—b and has two apses. Let 
ab = k' 2 , and choose arbitrarily k>a>a, and P=k 2 l a . Then the orbit lies 
between w=« and u—P, and has Property II. 

Case II. When//M 3 always increases with u, select k and « arbitrar- 
ily [<*>&], and let P=k 2 /a. Then the vanishing of du/d o for u= a , P, deter- 
mines a value of h, h 2 =P(u, )/u, , such that P(*)/a<h l <P(P)/P. Thus the 
trajectory has two branches of the kind considered in Theorem II. By 
varying « a family is obtained having Property II. 

8. Admitting Property III. Not all forces admit orbits in which 
the angle from the radius vector of a given length to the nearest pericenter 
is V2. For example, the laws represented by f—kr n where n>l have their 
apsidal angles all less than */2.* But for any force, /, such that u.f is an 
increasing function of u, all apsidal angles are greater than */2; and by 
Corollary 1, if also f/u 3 is everywhere a decreasing function, any pair of 
apsidal values « and P are admissible. 

Let c be an arbitrarily chosen constant; and let 6 denote the angle 
from the radius vector where u=c to the nearest pericentral line, where 
u--=P>c. Then 



(15) o=f 



du 



i/{P 2 -u* + [F(u)-F(P) ] /h*Y 



Now for sufficiently great values of h, this integral differs arbitrarily little 
from 

'This and the following statement are established by a test contained in an unpublished paper by the writer, 
presented to the American Mathematical Society, December 30, 1908. 
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(16) f du/V^-u*), or V2-arcsin(c//J); 

hence there are values of h large enough to make 0<V2. But if « 
be chosen sufficiently near c [«<c], ©— <> may be made arbitrarily small [® 
denoting the apsidal angle]; so that, since ©>V2, h can be so chosen that 
o>k/2. Between this value of h and the "sufficiently great values," there 
must exist a value for which o—n/2. 

By varying P a family of orbits is obtained having Property III. For 
any value of P the necessary value of h is given by (15) on placing 0—n/2; 
and « is given by (11) when P and h are known. 

There is likewise a family having Property III with respect to the 
apocenters; for with a given «, h can be selected so small, that o shall be ar- 
bitrarily small; or P selected so near c that o shall be greater than V2. 
Hence, for a certain value of h, 0=--* /2. 

The determination of h in each case is unique; for in (15) if h be as- 
signed two different values h, and h t , hi >h lt the integrand involving h 2 is 
smaller throughout than that involving h L , so that only for one of these val- 
ues can 0=x/2. 

9. Admitting Properties I and II. The preceding sections show 
that no one of the three properties serves to characterize any particular law 
of force. But it is quite otherwise if laws be sought admitting both Prop- 
erties I and II. 

For if the constant of areas be proportional to the curvature at each 
apse; then h 3 —xp(P)=/*P(a), where x and > are constants. And if the 
family is to have aP=k 1 , x and c- must be equal, since as « and P each 
approach the value k, h s approaches xp(k) or ^P(fc). Therefore, also 
P(k*/fi)=P(fi); and from (11) follows 

(17) F(P)-F(k*/P) = (P i -k*/P*) [x.P(P)] % , 

which must hold as P varies through the range of values taken in the fami- 
ly. Differentiation of (17) with respect to P gives on replacing P(k*/P) by 
TO: 

(18) {l+k*/P*)P(P)=(P+k*/P 3 ) [X P0»)]* 

+ (V3) (P* -h*/P) [x P{P)]-*P'{P), 

or substituting kx for P and <l>(x) for x.p(kx) : 

uy > daf~V-ar Xk'x*-V 

a differential equation which is of Bernoulli's type, and is reducible to a lin- 
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ear equation by the substitution <t>=y~ 3 . A particular integral is y=(x* +1) 
/2*kx, while the solution of the auxiliary equation is y—cy/ix^ — 1)/2^ kx; 
so that the general solution is 

(20) ^-^=y=[ci/(x 4 -l)+x'+l]/2^kx, 
where c is an arbitrary constant; or 

(21) P(P)=8k e * 2 P a [c l / (P*-k*) +P* +k s ]~ 3 . 
Hence the most general law of force, f=u*P(u), would be 



(22) f=8k^" u 



[c|/0* 4 -& 4 )+tt 2 +fc 8 ] 8 ' 



But, whether c is real or imaginary [c^O], (22) gives imaginary values for 
/either when u>k or else when u<k; and such a force would not admit real 
orbits in both parts of the plane. Thus the only admissible value of c is 
zero; and (22) then reduces to the same form as (4). 

Hence laws of Type I are the most general which admit families of or- 
bits possessing Properties I and II; and these two properties serve, therefore, 
to characterize laws admitting circular orbits whose centers are elsewhere 
than at the center of force. 

§ 4. The Families for Particular Laws. 

10. For Newton's Law. Every orbit described under a force oper- 
ating according to Newton's law is a conic having a facus at the center of 
force. Expressed in terms of the constants of the conic, the force is given 
by 
(23) f=h 1 u i /a(l-e i ), 

so that, for all the orbits, h 3 /a(l-e 2 ) must be constant. If also h varies as 
the curvature at any apse, l/a(l-e 2 ) or a/6 8 , it follows that a(l— e 2 ), and 
hence also h, must be constant in the family. Hence, for the Newtonian law, 
the family having Property I is the totality of conies having a certain curva- 
ture at a vertex, or what is equivalent, having latera recta of a certain 
length. Evidently there is an infinitude of such families, obtained by 
changing the arbitrary constant. And further, each of these families pos- 
sesses Property III also. 

Again, to have the product of the apsidal distances constant requires 
the constancy of a 2 (1— e 8 ), or of the minor axis. Hence for Newton's law, 
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the family possessing Property II is the set of conies having minor axes of 
an arbitrary constant length. 

11. For the Law op the Direct Distance. All orbits described 
under a force varying directly as the distance, are conies whose centers lie 
at the center of force. The force, expressed in terms of the constants of 
the conies, is given by 

(24) f=*h*/a*b 2 u; 

so that h^/a'b 2 must be constant for all the orbits. If also h is proportional 
to a/b 2 , it is necessary that 6 6 , and hence also b, remain invariant. Thus, 
for this law, the family possessing Property I is the set of conies having 
minor axes of a given length. The constant of areas must vary as the maj- 
or axis in this family. To have Property I at the pericenters, interchange 
b and a in the preceding statements. Each of these families has Property 
III also. 

Finally, to have the product of the apsidal distances constant requires 
the constancy of a. b and also of h. Hence, for the law of the direct dis- 
tance, the family possessing Property II is the set of conies having the 
product of the axes constant, or what is equivalent for ellipses, having a 
given (arbitrary) area. 



ON THE IRREDUCIBILITY OF CERTAIN POLYNOMIALS. 



By JACOB WESTLUND, Purdue University. 



The object of the following note is to determine whether the two 
polynomials 

f\(x)-=(x—ai)(x—a s ) ... {x-a n )— 1 and 
fs (x)=(x— a,) (x-a 2 ) ... (x— a„)+l, 

where a,, a s , ..., a» are distinct integers, are reducible or irreducible. 

Let us first consider/, (x). If/, (x) were reducible we would have 

where <£(#) is irreducible and of a lower degree than n. Then since 

<t>(a,i)4'(<ii)=—l, i=l, 2, ..., n, 
we must have 



